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ABSTRACT 

A theoretical analysis of pool film boiling based on a cellular model is presented. The time-averaged 
cell configuration is postulated to adjust itself to maximize the rate of heat transfer. Theoretical ex- 
pressions for the heat-transfer coefficient, average vapor dome (bubble) radius, and optimum cell diameter 
(wavelength) were in close agreement with experiments. As a generalization of this work and others, a 
tentative nonequilibrium thermodynamic hypothesis is proposed that nature maximizes the entropy production 
in all situations characterized by macroscopic chaos. 

AUSZUG 

o 

g Eine theoretische Analyse einer siedenden Oberflache auf der Grundlage eines zellenartigen Modells 

fo wird gegeben. Die zeitgemittelte Zellenkonfiguration stellt sich voraussetzungsgemass so ein dass sie die 
w Geschwindigkeit der Warmeubertragung maximiert. Ein theoretischer Ausdruck fur die Warmeleitzahl, den 
Darapfdomradius (Blase) und den optimalen Zellendurchmesser stimmte sehr nahe mit Versuchen uberein. Als 
eine Verallgeraeinerung dieser Arbeit fluid Anderer wird eine vorlaufige ausgleichslose thermodynamische 
Hypothese vorgebracht, namlich dass die Natur die Erzeugung der Entropie in alien durch makroskopisches 
Chaos charakterisierten Umstanden maximiert. 


AHHOTAUHH 

JJaH TeopeTHuecKHft aHajiH3 KHneHKH ruieHKH Ka ropH3 0HTajibHoft noBepxHOCTH b npeanojioaceHHH 
iiejiJiiojiHpHog MOflejm. IIpe^nojiaraeTCH uto ycpe^HeHHan bo BpeMeHH KOH(?)iirypaijHA srceftKH caMO- 
HacTpaHBaeTCH Ha MaKCHMajibHyio CKOpocTB Teiuionepe,uami. TeopeTHvecKHe BbipajKeHHH ajih koo(|xJ)h- 
mieHTa Tenjionepe^a^H, cpe^Hero pa^nyca ny3bipn h onTHMajibHoro ^naweTpa ( .hjihhu bojlhu ) 
H^eMKH naxo^HTCH b xoponieM corxacoBaHHH c oicciiepHMeHTajibHHMH ^aHHUMH . Ran o6o6meHHe stoh 
h 4pyrnx pa6oT ycxoBHo npe^jiaraeTCH rnnoTe3a TepMo^imaMHueeKOft HeypaBHOBemeHHOCTH npu 
KOTopoft npHpo^a caMa MHHHMH3HpyeT c03£aHHe shtpoiihh bo Bcex cjiyuanx MaKpocKonnuecKoro 
xaoca. 

INTRODUCTION 


Film-boiling heat transfer in pool systems is 
a subject of considerable interest, especially in 
the field of cryogenics, since relatively low vail 
temperatures will induce film boiling. A layer of 
dense fluid supported by a lighter one, for film 
boiling of liquid above vapor, is known to be an 
inherently unstable configuration in a gravita- 
tional field. At certain locations on the vapor- 
liquid interface, vapor will break through and es- 
cape under the influence of gravity into the bulk 
liquid. Photographs of film boiling show that 
these escape points are dome shaped and look like 
large bubbles and the surface is covered by a 
cell-type structure. 

Figure 1 illustrates the two-dimensional 
model used in the analysis. The model consists of 
a thin film of vapor between the wall and the sup- 
ported liquid with the vapor domes distributed sym- 
metrically over the surface of the plate. This 
model is similar to that used by Berenson [l];^ 
however, the analysis used in this paper is dif- 
ferent from that used by Berenson in that no re- 
course is made to experimental results or to one- 
dimensional small amplitude stability theory. 

The present paper takes the point of view 
that heat transfer is the dominating or driving 
force for both the flow and the instabilities; that 
is, the time-averaged configuration of the liquid- 
vapor interface is determined by thermal energy 
considerations. 
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BASIC MODEL AND EQUATIONS 

The major portion of the heat transport to a 
unit cell (see Fig. l) Is assumed to occur across 
the thin vapor film. The vapor domes are consid- 
ered so thick that essentially no heat is conducted 
into them. Physically, their function is to act as 
hydrodynamic sinks into which the generated vapor 
is dumped. The assumed symmetrical distribution of 
these sinks implies that the velocity field in the 
thin film is radially symmetric. The model pro- 
posed assumes that there is some time-averaged or 
ensemble -averaged configuration of the system, 
where all velocity, pressure, and temperature 
fields are at steady state and that this statisti- 
cally idealized configuration represents the aver- 
age behavior of the actual system. 

The heat- transfer coefficient can be obtained 
by solving the momentum and thermal energy equa- 
tions for flow and heat transport in the thin annu- 
lar vapor film belonging to a single cell. In or- 
der to solve the momentum equations the inertia 
terms in the Navier-Stokes equation are considered 
negligible. This assumption can be justified on 
the grounds that the small thickness of the film 
prevents the buildup of sizeable velocities. Addi- 
tional support is given to this assumption by 
noting in [2 and 3], which consider the Leiden- 
frost film boiling of drops, that the difference 
between the solutions obtained by retaining the 
inertia terms and by dropping them is shown to be 
negligible. In addition, the physical properties 
are assumed constant (but evaluated at the film 


^Numbers in brackets denote references. 
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temperature); radiation and radial temperature 
gradients in the energy relations are assumed neg- 
ligible. 

For the assumptions noted, the governing dif- 
ferential equations in cylindrical coordinates are 
as follows: 

Momentum: 


. dP / d 2 u 1 du 

0 = 'Be r- + H : + - 5" 


3r 


3P 


3r 2 r 3r 


3 2 v 




3r 2 


Continuity: 


du u dw 
dr r dz 


- 

(1) 

r 2 3z 2 / 

i 3w , a 2 v\ 

r Sz 2 / 

(2) 

= 0 

(3) 


Energy: 


w 


aT 

dz 



(4) 


The attendant boundary conditions are as fol- 
lows: 

(a) At the wall, z = 0. The velocity of the 
vapor is zero: 


u(r,0) = w(r,0) = 0 


( 5 ) 


The temperature of the wall is at a specified con- 
stant value: 


T(r,0) = T v (6) 

(b) At the liquid- vapor interface, 
z = 5 

The radial velocity vanishes: 

u(r,6) = 0 (7) 

The normal velocity has some constant value depen- 
dent on both momentum and energy considerations, 
and the surface is at the saturation temperature: 

w(r,S) = w & ss constant 

T(r,6) = T s (8) 

Since the bulk liquid is assumed to be at the satu- 
ration temperature T s , no heat is conducted into 
the bulk liquid. The momentum and energy equations 
are coupled at the interface by a balance of latent 
heat release and heat conduction through the vapor 
film to the evaporating surface: 

■ *" 6 <9) 

A static force balance requires that the average 
pressure of the vapor over the annular interface be 
equal to the system pressure P s at the flat por- 
tion of the liquid- vapor interface, els shown in Fig- 
ure 1: 


2n 



P(r,8)r dr = nP^R 2 - R 2 ) 


( 10 ) 


One further assumption is made about the pressure 
field in the next section to make the problem de- 
terminate. 

(c) Also, from radial symmetry, 

u(Rq,z) = 0 (11) 


METHOD OF SOLUTION 

A solution flow chart of the entire analysis 
is presented in Figure 2. Hiis chart contains the 
overall logic of the analysis and depicts how each 
subproblem ties in with the final results. Refer- 
ence to this chart after reading each subproblem 
section will aid in understanding the analysis. 


Stream Function 


The two-dimensional axisymmetric equations of 
motion with inertia terms neglected, often referred 
to as the "creeping motion" equations, are handled 
most easily by use of the stream function ^ de- 
fined by 


1 d\|r 
u = - 
r oz 

(12) 

1 d\{r 
r dr 

(13) 

The continuity equation is identically satisfied by 
these definitions. The momentum equations become 

E 4 t = 0 

(14) 

where E 2 is the operator 


2 a 2 i a a 2 

= 3r 2 ' r Sr * z 2 

(15) 

Solving Equation (14) for this problem 
separation of variables, results in 

, ty 

V = (*0 + A l z + *2^ + a 3 z3 ) 

(16) 

Substituting \Jr into Equations (12) and (13) 
and applying the boundary conditions (Eqs. (5) , 
(7), and (8)) result in the following expressions 
for the velocity components: 

. - 

(17) 

„ , . a!( r . . 8l) 

(18) 


where 0 is a real number and, at this point, 
completely arbitrary. 


Pressure Function 


The velocity expressions (Eqs. (17) and (18)) 
can be substituted into the Navier-Stokes equations 
to obtain 


= 5H p 2 (2z - 8) 

dz g c 


( 19 ) 
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These expressions can he integrated, subject to 
the boundary condition 


p (V s > = p 8- Vpi-p)t + f < 21 > 

to yield as a final expression for the pressure 



+ ^ + 3B&! (z 2 . 6z) (22) 

R 1 S c 


A derivation of Equation (21) appears in [l]. 
It suffices to say that the pressure in the dome is 
slightly less than the pressure P s at the liquid- 
vapor interface. This difference is due to the 
difference in static head corrected for surface 
tension and curvature. 

In order to determine the constant p 2 , the 
pressure function (Eq. (22)) is inserted into the 
force balance (Eq. (9)), which when integrated 
gives 



B^(Pj-p) 


jg_ f(l- f 2 ) 
g c 2 


ofy(l -f 2 ) 
1 


1 - f 4 (1 - f 2 ) 2 In f 
8 4 2 


(33) 


where f = /Rq. A criterion for evaluating Rq 
and R^ will be postulated in the optimum config- 
uration section. 

One of the main results of this analysis is 
that there is a unique relation between the rate 
of evaporation necessary to support the liquid and 
the thickness of the gap. This relation is ob- 
tained by substituting z = 8 into Equation (l7): 



Substituting for v from Equation (17), expanding 
the exponential function, integrating each term, 
neglecting terms of order 8 8 and higher, and 
solving for the gap thickness gives 


8 = 


2k(T v - T g ) 


n l/4 


K( x + To c p4J 


(27) 


Heat-Transfer Coefficient 


Once the gap thickness is known, it is a rela- 
tively simple matter to compute the heat-transfer 
coefficient, which is defined by 


dT 


k dz 

z=0 

< T w - 



(28) 


where hQ^ is the local heat-transfer coefficient 
based on the annular area surrounding the dome. 

The heat-transfer coefficient of interest, however, 
is based on the total area of the plate. If heat 
transfer in the dome itself is neglected, then 
h ce ii> the heat- transfer coefficient based on the 
total area of a cell, is related to h arm by 

Veil = “arm = Van* 1 ' p2 > (29) 

Employing the first derivative of Equation (25) at 
z = 0, Equations (27), (28), and (29), and neglect- 
ing terms of 8 8 and higher results in the follow- 
ing heat-transfer coefficient: 


w & = - | p 2 5 3 (24) 


Energy Equation. 

Integrating Equation (4) subject to the bound- 
ary conditions of fixed temperatures Tv and Tg 
at both the wall and the interface, respectively, 
yields the general solution (as differentiation 
will verify) 



where x and y are dummy variables. 

The interface energy balance is obtained by 
combining Equations (9), (24), and the first deriv- 
ative of Equation (25) evaluated at z = 8: 



This expression for the modified latent heat of 
vaporization was found by repeated application of 
the binomial theorem, while assuming that Cp AT/A 
was small. 


OPTIMUM CONFIGURATION OF LIQUID- VAPOR INTERFACE 

It will be shown in this section that there 
is an upper limit to the rate of heat transfer to 
a unit cell; that is, there exist "optimum" values 
(from the point of view of maximum heat transport) 
of the radii Rq and R-^, because very large 
cells and very small cells are incapable of maxi- 
mizing the heat transport. The heat transport to 
the entire system will be maximized if 
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^ h cell 

d h cell 

^ h cell 


^R-l 


k> 

df | 

% 

and 

^cell^^l^ 

Sh ceu/ R 0’ f ) 

se 3 

SRo 

R " ** 

R 1 

f = 4j 


= 0 


(32) 


(33) 


(These expressions were checked to insure that they 
did not lead to a minimum or a saddle point. ) Con- 
sider Equation (33) first. Taking the derivative 
of 0 2 , given by Equation (23), with respect to Bq, 
setting this equal to zero, and solving for Bq 
yield 
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gg c 

g(P z - p) 



(34) 


By substituting the value of Bq (Eq. (34)) 
into the expression for J3 2 (Eq. (23)) the heat- 
transfer coefficient given by Equation (30) becomes 


h cell “ 


k 3 X*pg(p z - p) 
“p(T w -\)l 


1/4 


0 


2 


0. 4716 


[(p z - p)g] 5/2 

p(ag c ) V2 


(41) 


In the top view of the cellular model (Fig- 
1(a)), small curved triangular regions exist on the 
plate that were not compensated for in the analy- 
sis. These curved triangular patches must be taken 
into account in order to evaluate the total heat- 
transfer rate. Since there are six triangular 
patches in contact with a single cell and three 
cells bordering each patch, the additional area 
attributed to a single cell is twice the area of a 
triangular patch. Drawing an equilateral triangle 
through the centers of three adjacent cells and 
evaluating the areas of the enclosed circular seg- 
ments readily shows that the area of the triangular 
patch is given by 


(42) 


The assumption is made that the local heat-transfer 
coefficient in the triangular patch region is the 
same as the heat -transfer coefficient in the annu- 
lar region surrounding a vapor dome. This assump- 
tion is based on the fact that the film thickness 
under the triangular patch will be equal to the 
film thickness under the annular area. Therefore, 


*tot = h ann(Vn + 2A A^ = h tot< A cell + 2A A^ (43) 


-11/4 


^ 7 i 


f 2 (i - f 2 ) 


3(6) (1 - f 2 ) 1 - f 4 . In f 

4 8 2 


(1-f 2 ) (35) 


where 


gc<* 


l 1 / 2 


[g(p z - p) 


(36) 


The derivative of h cell with respect to f, 
which was indicated by Equation (32), gives a com- 
plicated transcendental equation for the optimum 
value of f. It is much simpler to obtain the max- 
imum value of hceip numerically. The optimum 
value was found numerically to be 


Solving Equation (43) for h^^ and combining the 
result with Equations (29) and (42) yield 



The area ratios are evaluated from Equations (37) 
and (42); consequently, the heat-transfer coef- 
ficient (Eq. (38)) corrected for the triangular 
patches becomes 


htot = 0.410 


k 3 X pg(p z - p) 


P(T V - T S )I 


1/4 


(45) 


COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


f 


opt 


0.673 


(37) 


Thus, if Equation (37) is substituted into Equa- 
tion (35), the maximum value of the heat-transfer 
coefficient consistent with the proposed model is 


h ce n = 0.381 


fk 3 X*pg(p, - p)~j ^ 


p(T w - T g ) l 


(38) 


The optimum cell radius from Equation (34) is 
Rq = 3. 64 l (39) 

and the optimum dome radius from Equation (37) iB 
= 2.4 5 1 (40) 

while the value of p 2 is 


Equation (45) is identical in form to an ex- 
perimentally confirmed expression proposed by Ber- 
enson [l]. Kie only differences between the two 
results are the value of the prefactor constant 
and the definition of the modified latent heat of 
vaporization, but the lower value of the prefactor 
constant in the present results is compensated for 
by a higher modified latent heat, as shown in the 
following table: 


Parameter 

Berenson [l] ^ 

Present results 

Prefactor constant 

0.425 

0. 410 

Modified latent 
heat of vapori- 
zation 

4 + °- 5 ^) 

4*i¥) 


In Berenson's work, an empirical expression 
was used for the diameter of the dome. From a 
series of high-speed motion pictures of film 
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boiling of n- pentane and carbon tetrachloride, the 
diameters of the vapor domes that had just departed 
or were just about to depart were measured at var- 
ious heat fluxes. This characteristic dome diam- 
eter was correlated by the equation 


± 10 percent (46) 


This value is very close to the optimum vapor dome 
diameter given by Equation (40); namely. 


D 1 = 4.90 


g c q 

g(p} - p) 


DISCUSSION OF RESULTS 

Ofre formulation of the film boiling process in 
terms of a maximization principle used herein is a 
specific instance of what is called a "generalized 
optimal boundary value process. " Krzyvoblocki in a 
recent work [ 4 ] has formulated a very general 
statement of this kind of problem. In other words, 
all boundary value problems can be looked upon as 
control problems. The boundary conditions are the 
controls. The optimal boundary value problem is to 
select from all admissible boundary conditions a 
set for which some functional assumes an extreme 
(maximum or minimum) value. It is to be emphasized 
that the shape of the boundary is a boundary condi- 
tion. Apparently, certain kinds of physical situa- 
tions are naturally optimal boundary value proces- 
ses. Krzyvoblocki ( 4 ] gives some examples of this 
natural class of optimal or extremum processes, one 
of them being the dynamics of a liquid enclosed in 
a partly filled tank under conditions of zero 
gravity. Here, the unknown boundary conditions 
(distribution of vapor and liquid) can be deter- 
mined statistically by a variational principle. 

In the film boiling problem, since the vapor- 
liquid boundary is flexible, it is conceptually 
possible for the vapor liquid boundary to adjust 
itself to maximize heat transfer. In terms of the 
foregoing nomenclature, a generalized postulate on 
film boiling is proposed: Film boiling on a hori- 

zontal surface is an optimal boundary value process 
where the functional to be maximized is the total 
rate of heat transfer. Of all the possible con- 
figurations of the vapor-liquid interface, that one 
which maximizes the energy transport is realized. 

The final justification for the postulate is 
that results employing it are consistent with ex- 
perimental data, as shown in the preceding section 
of this report; however, the applicability of a 
maximum principle may even transcend this particu- 
lar application to film boiling. In [5] for ex- 
ample, Malkus advanced a novel theory on the nature 
of turbulent free convection. He proposed that in 
thermal turbulence the flow adjusts itself in such 
a way as to transfer the maximum amount of heat 
compatible with the boundary conditions. Without 
introducing experimental constants, a reasonably 
accurate prediction of the mean velocity distribu- 
tion for the turbulent flow was achieved. Other 
papers by Malkus [6 and 7] which treat finite- 
amplitude cellular convection and a theory of shear 
turbulence, respectively, also employ a maximiza- 
tion postulate. In addition, Zuber [8] has taken 
the Malkus theory of turbulent free convection as a 
basis for a general model of nucleate boiling. 


In the analysis it was tactily assumed that 
film boiling was fully established. The follow- 
ing question may arise: since nucleate boiling is 

a much more efficient process than film boiling and 
if nature is attempting to maximize the heat trans- 
fer, why is the heat -transfer rate degraded at some 
critical temperature difference (burnout) ? Appar- 
ently, these critical transition points must be 
looked upon as constraints. The variations of 
physical properties with temperature and pressure 
are thermodynamic constraints that always exist. 

Two examples illustrating how thermodynamic con- 
straints determine what heat -transfer regime will 
exist are as follows: 

(a) Above a certain wall temperature a 
liquid may be unable to wet the wall, which 
makes nucleate boiling impossible. 

(b) Above the thermodynamic critical pres- 
sure of a fluid, no boiling can occur. 

If Zuber' s theory [8] of boiling with the use 
of the Malkus model of turbulent free convection is 
correct, then a general statement on boiling pro- 
cesses can be formulated: Boiling heat transfer, 

both nucleate and film, is an attempt by nature to 
maximize the rate of heat transfer into the system, 
subject to the general constraints imposed by ther- 
modynamics and hydrodynamics. 

Which regime of boiling, nucleate or film, 
actually obtains in a given situation cannot be 
predicted from such a general postulate, but must 
be regarded as calculable from the general con- 
straints which are operating. However, given a 
regime and a choice among many alternatives (tur- 
bulence, cell distribution, etc.), nature will home 
in on that state of the system which maximizes the 
energy transport. 

On a philosophical basis, why is nature so 
obliging in maximizing the rate of heat transfer? 
The explanation lies in entropy considerations. In 
essence, everything said about maximizing the 
energy transfer can be recast in terms of an en- 
tropy production; that is, if the rate of heat 
transfer is being maximized, the rate of entropy 
production is also being maximized. If the hot 
wall is considered a heat reservoir at temperature 
T w and the bulk fluid is considered a heat sink at 
temperature T s , the rate of entropy production of 
the universe is then 

dS _ _Q_ 

dt - t b t w 

where Q is the total rate of heat transfer from 
the wall. Clearly, for fixed T w and T s , the 
rate of entropy production of the universe will be 
a maximum if the heat transport is a maximum. 

The entropy formulation is more basic and will 
probably comprehend more physical situations than 
the energy formulation. The generalization that 
suggests itself, then, is that in all chaotic-type 
situations the average behavior of the system is 
such that the entropy production is maximized. 

This generalization is the macroscopic analog of 
the Second Law of Thermodynamics. Clearly, the im- 
portance of this nonequilibrium thermodynamic pos- 
tulate is obvious, and its general validity war- 
rants further study and investigation. 


CONCLUSIONS 

A theoretical analysis of pool film boiling 
from a horizontal surface yielded the following: 
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1. A theoretical expression for the heat- 
transfer coefficient that agrees within 

4 percent of an experimentally confirmed ex- 
pression proposed by Berenson 

2. An optimum average vapor dome (bubble) 
radius that agrees within 5 percent of ex- 
perimentally observed breakoff diameters 

3- An optimum cell diameter (wavelength) which 
is intermediate between the critical and the 
most dangerous wavelength predicted by one- 
dimensional small -amplitude hydrodynamic 
stability theory 

In addition, a tentative nonequi librium ther- 
modynamic hypothesis was proposed which states that 
nature maximizes the rate of entropy production in 
all systems characterized by macroscopic chaos, 
such as turbulence, nucleate boiling, and film 
boiling. 


NOMENCLATURE 


Ad^ a i> 

Ag,A 3 

^ann 

A cell 

A A 

CpjCg 

C P 

% 

% 

E 2 

f 

B 

®c 

h ann 

hcell 

h tot 

k 


s 

Q 

^tot 

% 

% 

r 

S 

T 


w 

u 

w 

w 6 

x 

y 


arbitrary constants 

areas of annular cell region, sq cm 
area of cell, sq cm 

area of curved triangular patch, sq cm 
arbitrary constants 

vapor specific heat at constant pressure, 
cal/(g)(°K) 

breakoff diameter of vapor dome, cm 
diameter of cell, cm 
diameter of vapor dome, cm 

differential operator defined by Eq. (15) 
ratio of vapor dome to cell radius, R±/R$ 

local value of gravity, cm/sec 2 
conversion factor in Newton* s law of 

motion, 1.0 g-cm/( sec 2 ) (dyne) 
heat- transfer coefficient based on annu- 
lar area, cal/(sq cm)(sec)(°K) 
heat- transfer coefficient based on cell 
area, cal/(sq cm)(aec)(°K ) 
heat-transfer coefficient based on total 
plate area, cal/(sq cm) (sec) (°K) 
vapor thermal conductivity, cal/(cm) (sec) 

(°K) 

characteristic length given by Eq. (36), 
cm 

absolute pressure, dyne/sq cm 
pressure of system, dyne/sq cm 

rate of heat transfer, cal/sec 
total rate of heat transfer from plate, 
cal/sec 

radius of cell, cm 
radius of vapor dome, cm 
radial coordinate, cm 
entropy, cal/°K 
temperature, °K 
T w - T s , °K 

saturation temperature, °K 
wall temperature, °K 

radial velocity, cm/sec 
axial velocity, cm/sec 
evaporation velocity, cra/sec 

dummy variable, cm 
dummy variable, cm 
axial coordinate, cm 


a thermal diffusivity, k/pCp, cm/sec 2 
0 2 constant given by Eq. (41), l/(sq cm) (sec) 

5 thickness of vapor film, cm 

A latent heat of vaporization, cal/g 

A* modified latent heat of vaporization, cal/g 
M- vapor viscosity, g/(cm)(sec) 

p vapor density, g/cu cm 

Pi liquid density, g/cu cm 

a surface tension, dyne/cm 

i stream function, cu cm/sec 


REFERENCES 


1. Berenson, P. J. : "Film-Boiling Heat Transfer 

from a Horizontal Surface," J. Heat Transfer, 
vol. 83, no. 3, August, 1961, pp. 351-358. 

2. Baumeister, K. J. ; T. D. Hamill; F. L. Schwartz; 

and G. J. Schoessow: "Film- Boiling Heat 

Transfer to Water Drops on a Flat Plate," 

AIChE Preprint 24, presented at the ASME- 
AIChE Eight National Heat Transfer Conference, 
Los Angeles, August 8-11, 1965. 

3. Baumeister, K. J. ; and T. D. Hamill: "Creep- 

ing Flow Solution of the Leidenfrost Phenome- 
non," NASA TN D-3133, December, 1965. 

4. Krzywoblocki , M. Z. v. : "Certain Epihydrody- 

namic Problems as Generalized Optimal Boundary 
Value Processes," ASME Paper 65-AV-36, pre- 
sented at the ASME Aviation and Space Confer- 
ence, Los Angeles, March 14-18, 1965. 

5. Malkus, W. V. R. : "The Heat Transport and Spec- 

trum of Thermal Turbulence, " Royal Soc. Pro c. , 
Ser. A, vol. 225, no. 1161, August 31, 1954, 
pp. 196-212. 

6. Malkus, W. V. R. ; and G. Veronis: "Finite 

Amplitude Cellular Convection, " J. Fluid 
Mech. , vol. 4, pt. 3, July, 1958, pp. 225- 
260. 

7. Malkus, W. V. R. : "Outline of a Theory of Tur- 

bulent Shear Flow," J. Fluid Mech. , vol. 1, 
pt. 5, November, 1956, pp. 521-539. 

8. Zuber, N. : "Nucleate Boiling. The Region of 

Isolated Bubbles and the Similarity with 
Natural Convection," Intern. J. Heat Mass 
Transfer, vol. 6, January, 1963, pp. 53-78. 




(a) Top view: Showing symmetrical distribution of cells 
on surface. 



Side view: Showing direction of vapor flow into domes, 
axial coordinate z, radial coordinate r, cell radius Rq, 
dome radius Rj, and gap thickness 6. 

Figure 1. - Idealized cellular model for film boiling on horizontal surface. 









